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Abstract 

We explore cosmology in the decoupling limit of a non-linear covariant ex- 
tension of Fierz-Pauli massive gravity obtained recently in arXiv:1007.0443. 
In this limit the theory is a scalar-tensor model of a unique form defined by 
symmetries. We find that it admits a self-accelerated solution, with the Hub- 
ble parameter set by the graviton mass. The negative pressure causing the 
acceleration is due to a condensate of the helicity-0 component of the massive 
graviton, and the background evolution, in the approximation used, is indis- 
tinguishable from the ACDM model. Fluctuations about the self-accelerated 
background are stable for a certain range of parameters involved. Most sur- 
prisingly, the fluctuation of the helicity-0 field above its background decouples 
from an arbitrary source in the linearized theory. 

We also show how massive gravity can remarkably screen an arbitrarily 
large cosmological constant in the decoupling limit, while evading issues with 
ghosts. The obtained static solution is stable against small perturbations, 
suggesting that the degravitation of the vacuum energy is possible in the 
full theory. Interestingly, however, this mechanism postpones the Vainshtein 
effect to shorter distance scales. Hence, fifth force measurements severely con- 
strain the value of the cosmological constant that can be neutralized, making 
this scheme phenomenologically not viable for solving the old cosmological 
constant problem. We briefly speculate on a possible way out of this issue. 



1 Introduction and summary 



The observed late-time acceleration of the Universe [1], and the Cosmological Con- 
stant problem (see reviews [2, 3]), remain two of the most tantalizing, mutually 
connected puzzles at the interface of particle physics and cosmology. 

A promising approach to the late-time acceleration enigma is to invoke new 
degrees of freedom, belonging to the gravitational field itself (as in massive gravity), 
that give rise to the cosmic speed-up. This framework postulates the existence of a 
new energy scale - set by the graviton mass - which is very low; nevertheless, this 
scale is technically natural in the quantum-field theoretical sense. This approach, 
as is known by now, is challenging theoretically (hence, is interesting) , and happens 
to have robust observational predictions. 

Such a scenario was first worked out in a context of the DGP model [4] in Refs. 
[5, 6], where the cosmic acceleration is due to the helicity-0 component of a five- 
dimensional graviton. Hence, the solution is said to be sclf-accclerating. 

Regretfully, in the context of DGP, the self-accelerating solution is plagued by 
negative energy ghost-hke states in the perturbative approach [7, 8, 9], and despite 
the issue of whether or not the negative energy perturbations could be continued 
in the full nonlinear theory [10], the existence of non-perturbative negative- energy 
solutions [7, 11, 12] makes the self-accelerating branch unsatisfactory (in spite of the 
interesting finding of Ref. [13] that the quasi-classical approach does not seem to 
reveal the instabihties of this solution) . 

Certain generahzations of the DGP model, however, allow for stable self-accele- 
rating solutions, either by constructing an explicit braneworld model [14] where the 
negative energy ghost disappears, or by extending the decoupling limit of DGP to 
the "galilean" invariant interactions, [15]. 

In this work, we show for the first time that a theory of massive gravity may 
produce a self- accelerated geometry while being free of the problems that arise in the 
self-accelerating branch of DGP. In particular, we will work in a certain approxima- 
tion in which the helicity ±2, ±1, and helicity-0 modes of the massive graviton de- 
couple from each other in the linearized theory, while the nonlinear self-interactions, 
and interactions between them, are captured by a few leading higher-dimensional 
terms in the Lagrangian; this approximation constitutes the decoupling limit. 

In this approximation, we will show the existence of the self-accelerated solution, 
around which small fluctuations are stable. The acceleration is due to a condensate 
of the helicity-0 fleld, which in the decoupling limit is reparametrization invariant. 
On the other hand, since the helicity-0 is not an arbitrary scalar, but descends from 
a full-fledged tensor fleld, it has no potential, but enters the Lagrangian via very 
speciflc derivative terms flxed by symmetries [16]. These terms generate the negative 
pressure density which causes the accelerated expansion with stable fluctuations, as 
will be discussed below. 

From the observational point of view, the obtained self-accelerating background 
is indistinguishable, in the approximation used, from that of the ACDM model. As 
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to the fluctuations, however, the hehcity-0 could have introduced some differences. 
For instance, at cosmological distance scales it could have given an additional force 
leading to, e.g., changes in the growth of structure [17, 18] , while at shorter scales still 
being strongly screened via the Vainshtein mechanism [19], guaranteeing the recovery 
of General Relativity with tiny departures [19, 20], which may also be measurable 
[21, 22] in high-precision Laser Ranging experiments [23] (for recent detailed studies 
of the Vainshtein mechanism see Refs. [24]). All the above takes place in the DGP 
model. However, this is not what happens on the self-accelerated background in 
the massive theory: Surprisingly-enough, the fluctuation of the helicity-0 on this 
background decouples in the linearized approximation from an arbitrary source! 
Thus, the astrophysical sources need not excite this fluctuation, in which case one 
recovers exactly the ACDM results. It is likely, however, that this similarity of 
the self- accelerated solution and its fluctuations to the ACDM results will not hold 
beyond the decoupling limit {i.e. will not hold for the horizon-size scales). 

Furthermore, if we wish to tackle the Cosmological Constant problem (CCP), S. 
Weinberg's no-go theorem makes it impossible to find dynamical solutions within 
General Relativity (GR) without involving fine-tuned parameters, [2]. The idea of 
infrared (IR) modification of gravity, however addresses this puzzle by accepting a 
large vacuum energy and modifying instead the gravitational sector in the IR, so that 
vacuum energy gravitates very weakly [25] . Such a source would not manifest itself 
as strongly as naively anticipated in GR, i.e. it would be degravitated, while all the 
astrophysical sources would exhibit the GR behavior [26]. As shown in Ref. [25, 26], 
one can think of degravitation as a promotion of Newton's constant to a high pass 
filter operator thereby modifying the effect of long wavelength sources such as a CC 
while recovering GR on shorter wavelengths. In particular, theories of massive and 
resonance gravitons exhibiting the high pass filter behavior to degravitate the CC 
[25]. Moreover, it was shown in Ref. [27] that any causal theory that can degravitate 
the CC is a theory of massive and resonance gravitons. 

It is important to emphasize that in theories of massive gravity degravitation is 
a causal process (unlike more general theories considered in [26]). The real measure 
of whether or not a source is degravitated is given by its time evolution. During 
infiation for instance, the vacuum energy driving the acceleration of the Universe will 
not be degravitated for a long time. It is only after long enough periods of time that 
the IR modification of gravity kicks in and can effectively slow down an accelerated 
expansion [25, 26]. Hence, a crucial ingredient for the degravitation mechanism to 
work is the existence of a (nearly) static solution in the presence of a cosmological 
constant towards which the geometry can relax at late time (or after some long 
period of time). Indeed, Ref. [27] studied linearized massive gravity demonstrating 
that in this approximation degravitation takes place after a long enough period of 
time. 

In this paper, we focus on the hard mass case using the generalized Fierz-Pauli 
theory of massive gravity, as derived in [16]. We show that this model allows for 
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static solutions while evading any ghost issues at least in the decoupling limit. In 
this framework an arbitrary vacuum energy can be neutralized by the effective stress- 
tensor of the helicity-0 component of the massive graviton. Small fluctuations around 
this solution are shown below to be stable, as long as this static solution exists. 

Moreover, we find that the energy scale at which the interactions of the hclicity- 
modes become nonlinear is affected by the scale of the dcgravitated cosmological 
constant - the interaction scale being higher for larger values of the CC^. 

On the one hand, it is intriguing that the interactions of the helicity-0 can be 
kept linear up to the energy scale which is significantly higher than what it would 
have been in a theory without the CC. However, this very same phenomenon also 
creates a problem by postponing Vainshtein's recovery of GR to shorter and shorter 
distance scales. As a result, the tests of gravity impose a stringent upper bound on 
the vacuum energy that can be degravitated in this framework without conflicting 
measurements of gravity. Disappointingly, this upper bound turns out to be of the 
order the critical energy density of the present-day Universe, (10~^ eV)'^ - the value 
that does not need to be degravitated. 

A possible way out of this difficulty may be to envisage a cosmological scenario 
in which degravitation of the vacuum energy takes place before the Universe enters 
the radiation dominated epoch - say during the inflationary period, or even earlier. 
By the end of that epoch then the cosmology should reset itself to continue evolution 
along the other branch of the solutions that exhibits the standard early behavior 
followed by the self-acceleration, found in the present work. The existence of such a 
transition would depend on properties of the degravitating solution in the full theory. 
Since we have no detailed knowledge of this solution at the time of this writing, we 
have no concrete mechanism to substantiate the above scenario. Therefore, in what 
follows we will not rely on it. Instead, we emphasize that there still are two important 
virtues of the degravitating solution with the low value of the degravitated CC: (I) It 
is a concrete example of how degravitation could work in four-dimensional theories 
of massive gravity without giving rise to ghost-like instabilities. (II) As we will 
show, the degravitated solution with small values of CC can be combined with the 
self-accelerated solution discussed above, to give a satisfactory solution that is in 
agreement with the existing cosmological and astrophysical data. 

Last but not least, the solutions found in the decoupling limit do not necessarily 
imply the existence of the solutions with identical properties in the full theory. Nev- 
ertheless, the decoupling limit solutions should capture the local dynamics at scales 
well within the present-day Hubble four- volume, as argued in [15]. On the other 
hand, at larger scales the full solutions may be very different from our ones. These 
differences would kick in at scales comparable to the graviton Compton wavelength. 
Therefore, our solutions should manifest themselves at least as transients lasting 
long cosmological times. 

Organization of the paper is as follows. In section 2, we review the generalized 

"'^In this work wc will use interchangeably the notions of vacuum energy and CC, although there 
could be a big difference between the two when it comes to IR modified gravity [25]. 
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Fierz-Pauli theory of massive gravity and discuss its ghostless decoupling limit. We 
then start by focusing on self- accelerating solutions in section 3, first deriving the 
background solutions, then testing their stability, and finally studying the implica- 
tions for late-time cosmology. We then explore the cosmology in the presence of 

a cosmological constant in section 4, proving the existence of a stable dcgravitat- 
ing branch of solutions, and analyzing the stability of the dc Sitter branch. Brief 
discussions of the degravitating solution are given at the end of section 4. 

2 The Formalism 

Search for a consistent theory of a massive spin-2 field goes back to the original 
work of Fierz and Pauli [28]. Whereas any massive gravity should reduce to the 
Fierz-Pauli (FP) theory at the quadratic level [29], a generic nonlinear extension 
exhibits the sixth degree of freedom - the so-called Boulwarc-Dcscr (BD) ghost [30]. 
This sixth mode produces severe instabilities on cosmological backgrounds [31], as 
well as on locally nontrivial asymptotically fiat backgrounds (such as that of a point 
source, for instance) [32, 33, 34]. 

This problem is usually related to the helicity-0 sector of massive theories [32]. 
The latter can efficiently be studied in the decoupling limit, where the sixth mode is 
hidden in higher-derivative nonhnear terms for the hehcity-0 [32, 33, 34]. Such terms 
make the Cauchy problem ill-defined, unless additional initial data are supplied. 
This corresponds to an additional, sixth, degree of freedom which shows up as a 
ghost-like linear mode on various backgrounds mentioned above. 

Up until recently it was thought that the cancellation of the higher-derivative 
nonlinear terms for the helicity-0 was not possible [33]. However, recently an explicit 
construction was given in Ref. [16] in which all the nonlinear terms for the helicity-0 
with more that two time derivatives cancel. Below we briefly review these results 
and recast them in a more convenient form. We refer to Ref. [16] for more detailed 
discussions. 

Consider a 4D covariant theory of a spin-2 field [32], which, once expanded on 
Minkowski space-time gives a graviton of mass m: 



L = Ml,^gR - ^^v^ {U2{g. H) + Us{g, H) + U^{g, H) + U^{g, H) ■■■).{!) 



Here UiS denote the mass and potential terms of i order in 



Us{g,H) = c,Hl^ + C2HHl + c^H\ 

U,{g,H) = d^Hi, + d2HHl, + dsHlHl^ + d,H'Hl, + d,H\ 



(2) 
(3) 
(4) 



(5) 
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Index contractions arc performed using the inverse metric g'^"; the coefficients Cj, di 
and fi are a priori arbitrary. The tensor H^^, is not an independent entity; it is 
related to the metric tensor as i?^,^ = g/^i, — rjahd^(p°'dy(p^ , where a,b — 0, 1, 2, 3, 
r)ab — diag(— 1, 1, 1, 1), and Hf^^, is a covariant tensor as long as the four fields 
transform as scalars under a change of coordinates [32]. Hence, the potential terms 
in (1) can be rewritten as functions of the metric g and the specific combination of 
the four scalars 93", as U{g, S), where S^,^ = rjabdnip'^di.ip'^. However, we will not be 
exploiting the latter representation in the present work. Instead, following [32] we 
expand in terms of the coordinates x", and the field tt", as = (x" — tt'*) 6^ , 
and using the convention, g^i, — rj/j^i, + h/j^i^/Mpi, we obtain 

= 77^ + ^M^i- + ^'^^M ~ Vai3d^,n"d^Tr^ . (6) 

The TTo's represent the Stiickelberg fields that transform under reparamctrization to 
guarantee that the tensor H in (6) transforms covariantly. In the unitary gauge one 
could put TTa — (or, 93" = x"S^), in which case (1) reduces to the standard FP 
theory extended by a potential for the field /i^j,. However, this is not a convenient 
way of dealing with these degrees of freedom. Instead, it is more instructive to retain 
tTq and fix a gauge for h^^. 

The theory (1) was studied in detail in [16, 35], and a two-parameter family of 
the coefficients was identified for which no sixth (ghost) degree of freedom arises 
in the decoupling limit^. In these theories the higher derivative nonlinear terms 
either cancel out, or organize themselves into total derivatives. For these ghostless 
theories, the decoupling limit is defined as follows^ 

m ^ 0, Mpi 00, A3 = {Mpim^y/^ fixed. (7) 

In what follows, we will focus on the hcIicity-2 and helicity-0 modes, and ignore the 
helicity-1 modes as they do not couple to a conserved stress-tensor at the linearized 
level, and, therefore, can be set to zero self-consistently (see, however, important 
comments on this at the end of section 3.2). 

We therefore use the following decomposition for H^j, in terms of the canonically 
normalized helicity-2 and helicity-0 fields after setting tTq = 9a7r/A| 

- m;; + ^^i Af — ■ 

Then, one can show by direct calculations [16] that the Lagrangian (1) reduces in 
the decoupling limit to the following expression 

>c = -W^s';i^Kp + h^"" E 73^4"^^"], (9) 

n=l -^3 



^Interestingly, a recently proposed extension of General Relativity by an extra auxiliary di- 
mension [36, 37], automatically generates the coefRcients from this family at least up to the cubic 
order. 

^By "ghostless" we mean a theory with no ghost at least in the decoupling limit, implying that 
even if the BD ghost exists in the full theory, it must have a mass larger than the scale A3, [16]. 
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where the first term represents the usual kinetic term for the graviton, Oi = —1/2, 
and 02,3 are two arbitrary constants, related to the two parameters from the set 
{ci,di} which characterize a given ghostless theory of massive gravity. The ex- 
pression {Eh) denotes the linearized Einstein operator acting on h^i, defined in the 



standard way: S'^^Kp = -\{nh^^-d^daK-d^dah'l^d^d^h-r}^,^nh^r]f,^dadph'^^). 



The three symmetric tensors X/t" [11] are composed of the second derivative of 
the helicity-0 field II^i^ = d^di^n. The vr field ends up being gauge invariant in 
the decoupling limit. In order to maintain reparametrization invariance of the full 
Lagrangian the tensors Xy^"'*[n] should be identically conserved. These properties 
uniquely determine the expressions for xju^ at each order of non-linearity. The 
obtained expressions agree with the results of the direct calculations of Ref. [16]. A 
convenient parametrization for the tensors xj^^ which we adopt in this work is as 
follows: 



The remarkable property of (9) is that it represents the exact Lagrangian (excluding 
the helicity-1 part) in the decoupling limit: All the higher than quartic terms vanish 
in this limit, making (9) a unique theory to which any nonlinear, ghostless extension 
of massive gravity should reduce in the decoupling limit [16]. 

If external sources are introduced, their stress-tensors then couple to the physical 
metric hf^iy. In the basis used in (9) there is no direct coupling of vr to the stress- 
tensors. Hence, the Lagrangian (9) is invariant w.r.t. the shifts, and the "galilcan" 
transformations in the internal space of the tt field, d^Tr d^n + v^, where is 
a constant four-vector. The latter invariance guarantees that there is no mass nor 
potential terms generated for vr by the loop corrections. 

The tree-level coupling of vr to the sources arises only after diagonalization: The 
quadratic mixing h'^'^xju) , and the cubic interaction h'^'^xju) , can be diagonalized 
by a nonlinear transformation of /i^,^, that generates the following coupling of tt [16] 



Moreover, the above transformation also generates all the Galileon terms for the 
helicity-0 field, introduced in a different context in Ref. [15]^. 

Since the Galileon terms are known to exhibit the Vainshtein recovery of OR 
at least for static sources [15], so does the above theory with as = 0. The quartic 

'^The usual coupling of the Galileon field to the stress-tensor, ttT , considered in the generic 
Galileon theories, violates the "galilean" invariance that is needed to protect the Galileon field 
from acquiring a mass and potential terms. Here however, the invariance is manifest before diag- 
onalization and the theory is protected. 



xO-J[U]^e^-^s/^,U^0, 

x(2)[n] = £/''%/%n„;3np. 



(10) 
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interaction h^'^XlJ , however, cannot be absorbed by any local redefinition of /i^j/. 
It is still expected though to admit the Vainshtein mechanism. 

However, as we will show in the next section, on the self- accelerated background 
the fluctuation of the helicity-0 field decouples from an arbitrary source, making the 
predictions of the theory consistent with GR already in the linearized approximation. 
This decoupling is a direct consequence of the self- accelerated background and the 
specific form of the coupling (11). 



3 The Self- Accelerated Solution 

The universality of the decouphng hmit Lagrangian (9) for the class of ghostless mas- 
sive gravities, suggests the possibility of a fairly model-independent phenomenology 
of the massive theories that should be captured by the limiting Lagrangian (9). In 
the present section, we will be interested in the cosmological solutions in these theo- 
ries. We will directly work in the decoupling limit, which implies scales much smaller 
than the Compton wavelength of the graviton. In the case of the self- accelerated de 
Sitter solution for instance, this corresponds to probing physics within the Hubble 
scale, which as one would expect, is set by the value of the graviton mass. 



3.1 The solution in the decoupling Umit 

Below we look for homogeneous and isotropic solutions of the equations of motion 
that follow from the Lagrangian (9). The helicity-0 equation of motion reads as 
follows: 

dc.d^h^'' {a,e,^f"^e/p, + 2^£/''-£/\H^ + S^s/'^^s/^^H^H.^) = , (12) 
while variation of the Lagrangian w.r.t. the helicity-2 field gives 

+ E tI^4"^ [n] = 0- (13) 

n=l -^3 

We are primarily interested in the self-accelerated solutions of the system (12)-(13). 
For an observer at the origin of the coordinate system, the de Sitter metric can 
locally (i.e., for times t, and physical distances much smaller than the Hubble 
scale H~^) be written as a small perturbation over Minkowski space-time [15] 

^ [1 _ lH-'x'^x^]r]^Ax^dxr (14) 

The linearized Einstein tensor for the (dimensionless) metric (14) is given by 

= -^r^'^"-^ = -Si^'r;... (15) 
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For the helicity-0 field we look for the solution of the following isotropic form 

TT = i qAlx'^Xa + bAlt + cAa , (16) 

where q, b and c are three dimensionless constants. 

The equations of motion for the hehcity-0 and helicity-2 fields (12)- (13), there- 
fore, can be recast in the following form 

(-l + 2a2q + 3a3qA =0, (17) 



2 

1 



MpiH^ = 2qKl 



-- + a2q + a^q^ 



(18) 



Solving the quadratic equation (17) for q (for H 0),we obtain the Hubble constant 
of the self-accelerated solution from (18). Its mag nitudc, _ A3/Mpi = m^, is 
set by the graviton mass, as expected (positivity of is one of the conditions that 
we will be demanding below). It is not hard to convince oneself that there exists 
a whole set of self- accelerated solutions, parametrized by 02 and a^. This range, 
however, will be restricted further by the requirement of stability of the solution, 
which is the focus of the next section. 

Before doing so, let us briefly analyze the four scalars (^". Using the ansatz, (16), 
their expression is given by 

(^" = (l-g)x"5«, (19) 

if we set 6 = 0. Thus the four scalars vanish in the special case of g = 1, and the 
metric is g^^i, = H^i,, so that the Lagrangian considered in (9) reduces to standard GR 
plus a CC (at least at the background level). This happens only if the parameters of 
our theory are such that 03 = ^ — |a2, which is not the regime we will be interested 
in - we will indeed show in what follow that the stability of the self-accelerated 
background implies q ^ 1- 



3.2 Small perturbations and stability 

Here we investigate the constraints that the requirement of stability imposes on a 
possible background. Let us adopt a particular solution of the system (17)-(18) and 
consider perturbations on the corresponding de Sitter background 

V = K'^ + Xnu, TT = tt'' + (/), (20) 

where the superscript h denotes the corresponding background values. The La- 
grangian for the perturbations (up to a total derivative) reads as follows 

C = -ix'^-£;f'Xa^ + 6(a2 + 3a3g)^^^^0n0-3a3^^^ 

+ ^^^X-X^[$] + |x-X(^[$] + ^, (21) 
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where $ denotes the four-by-four matrix with the elements <l>^j, = d^d^^cj). The first 
term in the first line of the above expression is the Einstein term for Xnvi the second 
term is a kinetic term for the scalar, and the third one is the cubic Gahleon. The 
second line contains cubic and quartic interactions between x^v and 0, which are 
identical in form to the corresponding terms in the decoupling limit on Minkowski 
space-time (9). None of these interactions therefore lead to ghost-like instabilities 
[16], as long as the (p kinetic term is positive definite. 

Most interestingly, however, there is no quadratic mixing term between x and (f) 
in (21). Since it is only x^v that couples to external sources T^i, in the quadratic 
approximation, then there will not be a quadratic coupling of (j) to the sources 
generated in the absence of the quadratic x ~ 4' niixing. Therefore, for arbitrary 
external sources, there exist consistent solutions for which the fluctuation of the 
hehcity-0 is not excited, (p — 0. On these solutions one exactly recovers the results 
of the linearized GR. The above phenomenon provides a mechanism of decoupling 
the helicity-0 mode from arbitrary external sources! This mechanism is a universal 
property of the self-accelerating solution in ghostless massive gravity. 

Hence, there are no instabilities in (21), as long as a2 + Sa^q > 0. The latter 
condition, along with the requirement of positivity of H^, and the equations of 
motion (17), requires that the following system be satisfied: 

1 2 
-- + 2a2q + 3a3q =0, 



MpiH^ = 2qAl 



2 1 
a2q + a3q - - 



> 0, 02 + Sa^q > 0, 



for the self- accelerating solution to be physically meaningful. The above system can 
be solved. The solution is given as follows 

2ao a, , , 

a2<0, -^<«3<-y, (22) 

while the Hubble constant and q are given by the following expressions 

= m^[2a2q' + 2a,e - > 0, ? = + ^^^^ . (23) 

3a3 3v2a3 

It is clear from (22), that the undiagonahzable interaction h'^'^X^J plays a crucial 
role for the stability of this class of solutions: All theories without this term [i.e. 
the ones with = 0) would have ghost-like instabilities on the self-accelerated 
background. Notice as well that in the regime (22) none of the scalars vanish, 
and our model therefore differs from GR with a CC. 

We therefore conclude that there exists a well-defined class of massive theories 
with the parameters satisfying the conditions (22), which propagate no ghosts on 
asymptotically flat backgrounds, and also admit stable self-accelerated solutions in 
the decoupling limit. 
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As we mentioned before, the helicity-1 field enters only quadratically, or in higher 
order terms in the Lagrangian, and hence, can consistently be set to be zero {i.e. it 
does not need to be excited by any other fields). Nevertheless, once a background 
configuration for the helicity-0 field is switched on, the higher-dimensional mixed 
terms of the helicity-0 and helicity-1 could in principle flip the sign of the Maxwell 
kinetic term, giving rise to a vector ghost that would only enter the Lagrangian 
quadratically or in higher powers; this field would couple to other fields at the 
nonlinear level. This certainly would not be a satisfactory state of affairs. 

By restoring back the helicity-1 field in our expressions, and performing direct 
calculations wc have shown that in the n^^ order in nonlinearities, where n < 6, the 
coefficient of the Maxwell term on the self-accelerated background is proportional to 
(— |-|-2a2g + 3a3g^), up the corrections that are of the (n-|-l)''^ order, [38]. Hence, up 
to these corrections, the Maxwell term vanishes on the self- accelerated background! 

If this were the full story we would get a theory of helicity-1 coupled infinitely 
strongly to the fluctuation of helicity-0 in the decoupling limit. However, quantum 
loop corrections will necessarily generate a nonzero Maxwell term, as it is not pro- 
tected by any symmetries. In these loops propagate the helicity-0 mode, as well as 
the matter fields to which the helicity-1 couples nonhnearly (for instance, one of the 
couplings being, duAad^AaT^""). 

Then, interpreting the value of the tree-level coefficient of the Maxwell term 
(which is zero) as an infinite value of the inverse of the running U{1) coupling at 
some UV scale Auv > A3, we obtain that at lower scales the coupling constant has 
a positive value as long as the theory is not asymptotically free (in other words, 
the U{1) couphng would have a Landau pole at some high scale Aj/y)^. Hence, the 
helicity-1 sector would not have a ghost, but the scale at which it would become 
nonhnearly interacting (the Vainshtein scale) would be parametrically (logarithmi- 
cally) smaller than A3. Since the helicity-1 field does not have to be excited by any 
source, this will not be a concern for us. 



3.3 Late-time cosmology 

In this subsection we discuss the relevance of the results obtained above for the 
late-time local cosmological evolution of the Universe. As seen from the decoupling 
limit Lagrangian (9), the helicity-0 mode tt provides an effective stress-tensor that 
is "felt" by the helicity-2 field: 



n=l ^3 



"2 + 02? + as? 



(24) 



It is this stress-tensor that provides the negative pressure density required to drive 
the acceleration of the Universe. Supplemented by the matter density contribu- 



^Alternatively, if the particle content is such that the theory has a negative beta function, then 
the infinite value of the coupling constant should be attributed to some far IR scale, A/^ <C A3, 
and at any scale greater than Aju the helicity-1 theory would have a finite positive coupling square. 
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tion, it leads to the usual ACDM - like cosmological expansion of the background 
in the sub-horizon approximation used here. This is clear form the fact that the 
stress-tensor (24) gives rise to a de Sitter background as was shown in the previous 
subsection. Hence, in the comoving coordinate system - which differs from the one 
used above - the invariant dc Sitter space will be the self- accelerating solution. 

All this can be reiterated by performing an explicit coordinate transformation 
to the comoving coordinates. This will be done in two steps. In the so-called Fermi 
normal coordinates, the FRW metric can be locally written in space and for all 
times, as a small perturbation over Minkowski space-time: 



ds^ ^ -[1 - {H + H^)x^]df + 



dx^ = (ry^, + hl^"^) dx^dx", (25) 



where the corrections to the above expression are suppressed by higher powers of 
H'^y^. The Fermi normal coordinates, on the other hand, are related to those used in 
(14) (in which the FRW metric is a small conformal deformation of Minkowski space- 
time), by an infinitesimal gauge transformation [15]. The latter does not change the 
expression (24), since T^^ is invariant under infinitesimal gauge transformations in 
the decoupling limit. On the other hand, the Fermi normal coordinates can be 
transformed into the standard comoving coordinates (tg, Xg) as follows [15] 



t - -H{t)x\ 



ait) 



1 + -H\t)x' 



(26) 



The stress-tensor of a perfect fluid, T^^^ = diag(p(tc), (i'^{tc)p{tc)5ij), transforms imder 
this change of coordinates (at the leading order in H'^y?) into the following expression 



p 

-H{p + p)x^ 



—H[p + p)x'' 



where all quantities in the latter expression are evaluated at time t. Note that the 
off-diagonal entries of the stress-tensor for the cosmological constant vanish in the 
Fermi normal coordinates, the same is true for T^j, as well. Hence, in all coordinate 
systems used the expressions for the stress-tensor on the self-accelerated solution is 
given by (24). 

Not surprisingly, the corresponding cosmological equations coincide with the 
conventional ones for the ACDM model, with the cosmological constant set by the 
mass of the graviton 



P 



a 
a 



(p+3p) + 



(27) 
(28) 



Here p and p denote the energy and pressure densities of matter and/or radiation, 
and = 6q [—^ + 02? + a^q^] is a constant that appears in (24). 
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As already mentioned, irrespective of the completion (beyond the Hubble scale) 
of the self-accelerated solution, it is locally indistinguishable from the ACDM model. 
At the horizon scales, however, it is likely that these two scenarios will depart from 
each other: As emphasized in the first section, the solutions found in the decoupling 
limit do not necessarily imply the existence of full solutions with identical proper- 
tics. Moreover, the decoupling limit Lagrangian (9) is derived from the full theory 
by dropping certain total derivative terms (see [16]), implying solutions that decay 
fast-enough at infinity. On the other hand, the solutions that we found in this sec- 
tion are given in the coordinate system where the fields grow at large distance/time 
scales. If these solutions are to be continued into the full theory, the latter should 
have an appropriate large scale behavior in this coordinate system. A given solution 
in the decoupling limit can just be a transient state of the full solution. Signifi- 
cant deviations of the latter from the former should kick in at distance/time scales 
comparable to the graviton Compton wavelength. 



4 Screening the Cosmological Constant 

4.1 Degravitation in generic theories of massive gravity 

One explicit realization of degravitation is expected to occur in massive gravity, 
where gravity is weaker in the IR, and the graviton mass could play the role of a high- 
pass filter [25]. In the approach of [25], the original theory was formulated as a higher 
dimensional model, the 4D reduction of which can be though as massive/resonance 
gravity with the mass term promoted into a specific differential operators determined 
by the underlying higher-dimensional construction. 

A more general approach was adopted in Ref . [27] , where the graviton mass was 
also promoted to an operator parameterized by a continuous parameter a 

m'(n) = mf ""^D", (29) 

the inverse graviton propagator is typically of the form 

^-1 ^ □ _ rnf^-'^^n^ , (30) 

so that for a < 1, gravity is weaker beyond wavelengths comparable to the graviton 
Compton wavelength mg ^. 

To take this mechanism a step further, the reliability of this argument within 
the non-linear regime is hence crucial. A trick to manifest the key interactions that 
arise in massive gravity is to work in the decoupling limit, where the usual GR 
interactions are suppressed, while the interactions of the new degrees of freedom 
are emphasized. This approach was first derived in Ref. [27], which we discuss first 
before turning to our considerations. As mentioned above, this limit is obtained 
by taking Mpi — )■ oo and m — > 0. However unlike in the decoupling limit of the 
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theory discussed in the previous section, the nonhnear dynamics in a generic model 
of massive gravity is governed by the scale 

^5-4a _ Mpim^(^-") . (31) 

In such models, it has been shown [27] that the helicity-0 (tt) and -2 (/i^i/) modes 
satisfy the following equations in the decoupling limit, 

- (3n(ni-%)^ + •••) = "mp^ ' ^^^^ 

where the physical metric is given by gi^i, = rj^u + (/'^i/ + 7r?7^,y)/Mpi. In the presence 
of a cosmological constant, = —^VtJ,i^j the solution for the helicity-2 mode is 

which is the usual GR solution. One can now check the condition for the existence of 
a (nearly) static solution towards which the geometry can relax at late times. In the 
language of the decoupling limit, this would happen if the helicity-0 mode compen- 
sates the helicity-2 mode contribution Trry^,^ = —h/^^ to maintain the geometry fiat 
Ofiu = Vfj-iy- However the configuration vr = Xx'^/6Mp\ is precisely the solution of (33) 
when the higher interactions vanish, i.e. GMpiDtt = — T = 8A. As shown in [27], 
such interactions cancel for tt ~ a;^ only if a < 1/2, hence implying that a generic 
theory of massive gravity amended with a nonzero CC can only have a static solu- 
tion when a < 1/2. In particular, in this language the DGP model [4] corresponds 
to a = 1/2 (see Ref. [7], but also [39]) hence explaining why this model does not 
bear static solutions with a brane tension, while promoting it to higher dimensions 
corresponds to a theory with a ^ for which the usual co dimension-two conical 
solutions can accommodate a tension without acceleration, [40, 41, 42, 43, 44]. 

The above results hold true for a generic theory of massive gravity. We now 
focus the analysis of the ghostless theory [16] reviewed in section 2, which strictly 
spealcing are not captured by the above a parametrization. The key difference in 
the ghostless case is that interactions for the helicity-0 mode are governed by the 
larger couphng scale A3 > A^. The form of these interactions in the ghostless theory, 
as well as the specific couplings to matter, play a crucial role in accommodating a 
degravitating branch of solutions, and this without being plagued by any instability 
at least in the decoupling limit. 

4.2 Degravitation in ghostless massive gravity 

For convenience we start by recalling the decoupling limit Lagrangian of (9) coupled 
to an external source 

^ = —h^'^^a, + h^"^ Yl T3?^4"^[n] + ^h^'^T,^ . (35) 

^ n=l ^3 ^^^Pl 
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The equations of motion for the hehcity-0 and 2 modes are then 



=1 -'^3 



and 



{a, + ^Dtt + ^ {[Ur - m) ) [ah - d^d.h-^ 



Baa 



(36) 



+^(a2n^.-3^(n 



2 



3a3 



(37) 



We now focus on a pure cosmological constant source, 7)^,^ 
of a similar ansatz as previously, 



TT 



2 A 

qx A 



2 A3 
3 ■ 



The equations of motion then simplify to 

n=l / 

(ai + 2a2g + Sa^q'^) = . 



A 



6Mi 



PI 



—Xrji^i,, and make use 

(38) 
(39) 

(40) 
(41) 



As we will see below, this system of equations admits two branches of solutions, a 
"degravitating" one, for which the geometry remains flat (mimicking the late-time 
part of the relaxation process), and a "de Sitter" branch which is closely related 
to the standard GR de Sitter solution. We start with the degravitating branch 
before exploring the more usual de Sitter solution and show that the stability of 
these branches depends on the free parameters 02,3, as well as the magnitude of the 
cosmological constant. 



4.2.1 The degravitating branch 

In this formalism, it is easy to check that the geometry can remain fiat i.e. H = 
and gi^j, = rin^, despite the presence of the cosmological constant. Such solutions are 
possible due to the presence of the extra hehcity-0 mode that "carries" the source 
instead of the usual metric. With H — 0, equation (41) is trivially satisfied, while 
the modified Einstein equation (40) determines the coefficient (which we denote by 
qo here) for the hehcity-0 field in (39), 
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in terms of the dimensionless quantity A = A/AgMpi. Notice that as long as the 
parameter as is present, Eq. (42) has always at least one real root. There is therefore 
a flat solution for arbitrarily large cosmological constant. 

Let us now briefly comment on the stability of the flat solution, as this has 
important consequences for the relaxation mechanism behind degravitation. We 
consider the field fiuctuations above the static solution, 

TT = ^goAg^^^'-^A, (43) 

Tf,y = -Ar^^^ + r^^, (44) 

where Qq is related to A via (42) and the coupling k is determined by 

K, = 2(ai + 2a2qo + SosQ'q) . (45) 
To the leading order, the action for these fluctuations is then simply given by 

£(2) ^ -^h^'^S^^^h^^ - ^h'^^XjlJ [$] + ^h^'^r,. , (46) 

with = d^dycj). The stability of this theory is better understood when working 
in the Einstein frame where the helicity-0 and -2 modes decouple. This is achieved 
by performing the change of variable, 

hf,iy = h^u + 4)7]^^ , (47) 

which brings the action to the following form 

£(2) ^ J_h^^S';:fKf, + ^-<t>n<j> + ^ [h^^ + <t>r^^^) r,. . (48) 

Stability of the static solution is therefore manifest for any region of the parameter 
space for which k is real and does not vanish. As already mentioned, if as 7^ there 
is always a real solution to (42), which is therefore stable for k 7^ 0. Furthermore, 
direct calculations to the 6**^ order show that the helicity-1 fluctuations will have a 
positive kinetic term as long as ^/{qo — 1) > 0. This suggests the presence of a flat 
late-time attractor solution for degravitation. The special case as = is discussed 
separately below. 



4.2.2 de Sitter branch 

In the presence of a cosmological constant, the held equations (40) and (41) also 
admit a second branch of solutions; these connect with the self-accelerating branch 
presented in section 3, and we refer to them as the de Sitter solutions. The param- 
eters for these solutions should satisfy 

oi + 2a2?dS + 3as?ds = > (49) 
A 2A^ 

^ds = ^T72- + TT ("i^ds + a2q% + as^ds) • (50) 
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This solution is closer to the usual GR de Sitter configuration and only exists if 
al > Sa^a^. The stability of this solution can be analyzed as previously by looking 
at fluctuations around this background configuration, 

TT = ^qdsAlx^ + (p, (51) 

hfMu = ~1^^6S.-^ Vnv Xnv 1 (52) 

T^^ = -Ar?^,. + T/,^ . (53) 
To second order in fiuctuations, the resulting action is then of the form 

£(2) ^ -l^^/^'^f^^^xa^ + ^^|p(a2 + 3a3to)0n</' + -^r^'r,^ ■ (54) 

It is interesting to point out again that the helicity-0 fiuctuation then decouples 
from matter sources at quadratic order (however the coupling reappears at the cubic 
order). Stability of this solution is therefore ensured if the parameters satisfy one of 
the following three constrains, (setting ai — —1/2 and A > 0) 

2al 1 - 3a2A - (1 - 2a2Xf^^ , , 

02 < and < < — — , (55 

3 - ^ 3A2 ^ ' 



or 



1 ^ 1 - 3a2A + (1 - 2a2A)3/2 
«2 < ^ and aa > , (56) 

or 

02 > ^ and as > (57) 

These are consistent with the results (22) found for the self-accelerating solution in 
the absence of a cosmological constant. Moreover, the requirement of stabihty of 
helicity-1 fiuctuations docs not impose further bounds on the parameters (see, dis- 
cussions at the end of section 3.2). Notice here that in the presence of a cosmological 
constant, the accelerating solution can be stable even when 03 = 0. This branch of 
solutions therefore connects with the usual de Sitter one of GR. 



4.2.3 Diagonalizable action 

In section 3 we have emphasized the importance of the contribution of xf} for 
the stability of the self- accelerating solution. However, in the presence of a nonzero 
cosmological constant, this contribution is not a priori essential for stability of either 
the degravitating or the de Sitter branches. Furthermore, since the helicity-0 and 
-2 modes can be diagonalized at the nonlinear level when aa = 0, as was explicitly 
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shown in [16], we will study this special case separately below. In particular, we will 
show that it leads to certain special bounds both in the degravitating and de Sitter 
branches of solution. 

Stability: To start with, when 03 — 0, the degravitating solution only exists if 

2a2~X < 3al . (58) 

This bound, along with the stability condition for the helicity-1 > 0, then 

also ensures the absence of ghost-like instabilities around the degravitating solution. 
Assuming that the parameters ai^2 = C>{1) take some natural values then the situa- 
tion a2 > implies a severe constraint on the value of the vacuum energy that can 
be degravitated. This is similar to the bound in the non-linear realization of massive 
gravity [37], as well as in codimension-two deficit angle solutions, A < m^M|j. The 
situation a2 < on the other hand allows for an arbitrarily large CC. 

On the other hand, the bound > 80103 for the existence of the de Sitter 
solution is always satisfied if 03 = 0. However, the constraints on the parameters 
(55) - (57) which guarantee the absence of ghosts on the de Sitter branch imply that 

2a2A > 3al . (59) 

In this specific case then, we infer that when the Sitter solution is stable, the de- 
gravitating branch does not exist, and when the degravitating branch exists the de 
Sitter solution is unstable. Therefore, at each point in the parameter space there 
is only one, out of these two solutions, that makes sense. In the more general case 
where 03 7^ the situation is however much more subtle and it might be possible to 
find parameters for which both branches exist and are stable simultaneously. 

Einstein's frame: Finally, to understand how this degravitating branch connects 
with the arguments in [27] and how it relates with Galilcon theories, let us now work 
instead in the Einstein frame, where the helicity-2 and -0 modes are diagonalized 
(which is possible as long as 03 = 0). The transition to Einstein's frame is performed 
by the change of variable [16, 45] 

V = h^u - 2ai7r77^^ + ^a^^Tr^^Tr , (60) 
such that the action takes the form 

jC = -V(^/.),. + 6a?7rnvr-^(a7r)^[n] + ^(a7r)^([nVm 

+ ^ V - 2ai7r77^, + ^a^Tra.TT^ (61) 

and the structure of the Galileon becomes manifest. Notice however, that the co- 
efficients of the different Galileon interactions are not arbitrary. Furthermore, the 
coupling to matter includes terms of the form df^7:d,^7:T'^^ , absent in the original 
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Galileon formalism [15]. Both of these distinctions play a crucial role in screening 
the cosmological constant - the task which was thought impossible in the original 
Gahleon theory. Here, however, as long as the bound (58) is satisfied, the solution 
for TT reads 

n^^qoA-lx'^ with Oi^o + 02^0 = (62) 
while the helicity-2 mode h^^ now takes the form 



with ^ being an arbitrary gauge freedom parameter. Fixing ^ = — 2a2?o-^3, 
physical metric is then manifestly flat: 

^ fj ^ . 2a2 



= riiiv - ( aigo + 02^0 + g I ^'^^lui^ + ]^ + 202^0^3) 

= 77m- (64) 

To reiterate, the specific nonlinear coupling to matter that naturally arises in the 
ghostless theory of massive gravity is essential for the screening mechanism to work. 
This allows us to understand why neither DGP nor an ordinary Galileon theory are 
capable of achieving degravitation. 



4.3 Phenomenology 

Let us now focus on the phenomenology of the degravitating solution. This mecha- 
nism relies crucially on the extra helicity-0 mode in the massive graviton. However 
tests of gravity severely constrain the presence of additional scalar degrees of free- 
dom. As is well known in theories of massive gravity, the helicity-0 mode can evade 
fifth force constrains in the vicinity of matter if the helicity-0 mode interactions are 
important enough to freeze out the field fiuctuations, [19]. 

Around the degravitating solution, the scale for helicity-0 interactions are no 
longer governed by the parameter A3, but rather by the scale determined by the 
cosmological constant A3 ~ (A/Mpi)^/'^. To see this, let us pursue the analysis of 
the fluctuations around the degravitating branch (43) and keep the higher order 
interactions. The resulting Lagrangian is then 
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with 



02 „a2 + 3a3go Mpi og 203 /MpiV 

F = -'^vp T' F = -Ap^l^J ' 

assuming 02,3 ~ To evade fifth force constrains within the solar system, the 

scale A should therefore be small enough to allow for the nonlinear interactions to 
dominate over the quadratic contribution and enable the Vainshtein mechanism. In 
the DGP model this typically imposes the constraint, A^/Mpi < (10~^^ eV)^, while 
this value can be pushed by a few orders of magnitude in the presence of Galileon 
interactions, [15, 46]. Therefore, the allowed value of vacuum energy that can be 
screened without being in conflict with observations is fairly low, of the order of 
(10-3 eV)^ or so. 

Notice that this maximal cosmological constant is at least of the same order of 
magnitude, if not better, than the tension that can be carried by a codimension-2 
brane embedded in six dimension with a Planck scale Mq. In this scenario, the 
maximal tension is of the order of A < 27rMg. From a four- dimensional point of 
view, this model with the brane-induced Einstein-Hilbert term looks like a theory 
of massive gravity with a graviton mass ~ M|/M|j. Phenomenology imposes 
the graviton mass to be m < lO^^^eV, which therefore imphes the upper bound of 
the brane tension, A < (lO^^ eV)^. 

The above constraint on the vacuum energy that can be degravitated makes the 
present framework not viable phenomenologically for solving the old cosmological 
constant problem. There may be a way out of this setback though: As mentioned in 
the first section, one may envisage a cosmological scenario in which the neutralization 
of vacuum energy takes place before the Universe enters the epoch for which the 
Vainshtein mechanism is absolutely necessary to suppress the helicity-0 fluctuations. 
Such an epoch should certainly be before the radiation domination. During that 
epoch, however, the cosmological evolution should reset itself -perhaps via some sort 
of phase transition - to continue subsequent evolution along the other branch of the 
solutions that exhibits the standard early behavior followed by the self-acceleration, 
found in the present work. This scheme would have to address the cosmological 
instabilities discussed in Refs. [47, 48]. Moreover, the viability of such a scenario 
would depend on properties of the degravitating solution in the full theory -which 
are not known. Therefore, we do not rely on this possibility. 

Nevertheless, there are certain important virtues to the degravitating solution 
with the low value of the degravitated CC. This is an example of high importance 
in understanding how S. Weinberg's no-go theorem can be evaded in principle. As 
already emphasized in [41, 42, 43, 44], such mechanisms evade the no-go theorem 
by employing a field which explicitly breaks Poincare invariance in its vacuum con- 
figuration TT ~ x^, while keeping the physics insensitive to this breaking. Indeed, 
physical observables are only sensitive to U^i, = d^di^i: which is clearly Poincare 
invariant, while the configuration of the tt field itself has no direct physical bearing. 
This is built in the specific Galileon symmetry of the theory, and is a consequence of 
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the fact that tt is not an arbitrary scalar field but rather descends as the helicity-0 
mode of the massive graviton. More precisely, under a Poincare transformation, 
x'^ — >■ A^x'' + a'*, the configuration for tt transforms as — >■ + t'^x'^ + c, with 
Vn — 2aiyA'j^ and c — which is precisely the Galileon transformation for tt under 
which the action is invariant. In other words the Poincare symmetry is still realized 
up to a Galilean transformation (or, there is a diagonal subgroup of Poincare and 
internal "galilean" groups that remains unbroken by the VEV of the tt field) . 

Thus, we have presented here the crucial steps towards a non-linear realization 
of degravitation within the context of massive gravity, and this, without introducing 
any ghosts (at least in the decoupling limit). The arguments presented here only 
rely on the decoupling limit and it is reasonable to doubt their validity beyond that 
regime. Fortunately, non-linear theories of massive gravity have been explicitly for- 
mulated in [36, 37], and static solutions in the fully non- linear regime have been 
presented in [37]. The absence of the ghost in theories of massive gravity requires 
the presence of additional symmetry projecting out the usual Boulware-Deser ghost, 
which can typically be thought as inherited from a higher dimensional fundamen- 
tal theory. It is therefore only natural to investigate massive gravity as arising in 
braneworld models embedded in (spurious) extra dimensions. The static solutions 
presented so far then embrace a much more physical meaning, where the quantity 
Hfj^i, plays the role of the extrinsic curvature on the brane, describing the brane posi- 
tion along the extra dimension(s). The fact that our model allows for flat solutions 
while carrying the cosmological constant with 11^,^ suggests that such models could 
be understood as flat branes embedded in extra dimensions, similarly as in [36, 37], 
[40, 41] and [42]. 

Some interesting work in Refs. [44] appeared during the completion of this 
manuscript. These have certain overlaps with the ideas of section 4 of the present 
work. In particular, Refs. [44] emphasize the role of Galileon flelds in the context of 
degravitation. These works differ, however, in several aspects from the present one. 
In particular, Ref. [44] rehes on the existence of two Gahleon flelds, as would arise in 
models with two extra dimensions, [41], whilst our model explores the degravitating 
solutions with a unique extra helicity-0 mode which naturally arises in the 4D theory 
of massive gravity. Our mechanism is possible thanks to the very speciflc coupling 
to matter that arises in a ghostless theory of massive gravity, and differ from the 
standard Galileon coupling. 
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